HOT SCALAR THEORY IN LARGE N: BOSE-EINSTEIN 
CONDENSATION * 



PETER ARNOLD 
reporting on work done in collaboration with 

BORIS TOMASIK 

Department of Physics, University of Virginia, P.O. Box 400714, Charlottesville, 

VA 22904-4714, USA 

I review the Bose-Einstein condensation phase transition of dilute gases of cold 
atoms, for particle theorists acquainted with methods of field theory at finite tem- 
perature. I then discuss how the dependence of the phase transition temperature 
on the interaction strength can be computed in the large N approximation. 



1 Phase Transitions in Hot Scalar Theories 

The standard example from particle physics of a scalar theory is the Higgs 
sector of electroweak theory, which has a phase transition (or a crossover) at 
a temperature of order the weak scale: say, a few hundred GeV or so. Let's 
focus on pure scalar theory by imagining setting the gauge coupling constant 
g w to zero. At finite temperature, the Higgs picks up a thermal contribution 
#AT 2 to its effective mass, and the effective potential becomes roughly of the 
form V{4>) ~ m 2 cS {T) 4> 2 + \<t> A with m 2 ff (T) = -M 2 + #AT 2 . At sufficiently 
high temperature, the #AT 2 turns the potential from being concave down at 
the origin to concave up, and so restores the symmetry that is spontaneously 
broken at zero temperature. 

Standard techniques for analyzing the phase transition between the hot, 
symmetry-restored phase and the cold, symmetry-broken phase are as follows, 
(i) Work in Euclidean time (for studying non-dynamical questions). The 
Euclidean time direction becomes periodic at finite temperature, with period 
j3 = 1/T. (ii) Near the second-order phase transition (T — > T c ) of the purely 
scalar model, the correlation length becomes infinite. Large distance physics 
becomes important, and there are large, non-perturbative, large-wavelength 
fluctuations, (iii) At large distances (Ek ~ k <§; T), the compact Euclidean 
time direction decouples, and one can reduce the original Euclidean theory to 
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a purely 3-dimensional effective theory of the zero-Euclidean-frequency modes: 



<Seff 



J d 3 x [\Vcf>\ 2 + m 2 cS \4>\ 2 + X cS \4>\ 4 + ■■■]. (1) 



(iv) Figure out what to do with the 3-dimensional theory (put it on a lattice, 
or whatever). 

2 Today's Talk: Bose-Einstein Condensation 

The purpose of today's talk is to show that the exact same techniques par- 
ticle physicists use to study relativistically hot scalar theories can also be 
used to study the Bose-Einstcin condensation phase transition of dilute gases 
of (for example) Rubidium atoms at T ~ 0.1 fiK. There's an identical 
three-dimensional effective theory to study the non-perturbative long-distance 
physics (Ek <C T) near the critical temperature: 



S cff = / d 3 x 



!■ 



V^ + r|^+ U |0| 4 + ... , (2) 



where I've switched to typical condensed-matter names (r and u) for the 
coefficients. For a non-relativistic problem, Ek ~ k 2 /(2m), and the long- 
distance condition Ek <C T for the validity of this effective theory becomes 
k < V2mT. 

The fact that three-dimensional effective theories can describe second- 
order phase transitions is old hat in condensed matter physics. What's dif- 
ferent about the dilute atomic gas problem is that, just as in weakly coupled 
relativistic problems, the coefficients of that effective theory are not phe- 
nomenological parameters but can be systematically matched to the micro- 
scopic physics! (This allows one to use the effective theory to calculate more 
than just universal properties of the transition.) 

3 Non-Interacting Non-Relativistic Bose Gas 

The path integral corresponding to a non-interacting, non-relativistic Bose 
gas in an external potential V(x) is 

Z = J[DiP}e l f dtL , (3) 
L = J d 3 x ip* ^id t + ^— V 2 - V(a;)^ tp + (chemical potential term). (4) 



where ip is a complex bosonic field and I've written the path integral, for the 
moment, in real time rather than Euclidean time. Why is this the path inte- 
gral? Note that the equation of motion, obtained by varying with respect to 
ip* , is just the Schrodinger equation (idt + V 2 /2m — V(x)^j ip = 0. The above 
path integral therefore describes the second quantization of the Schrodinger 
equation: it describes arbitrary numbers of particles, just like the standard 
path integral for QED describes arbitrary numbers of photons. In canonical 
quantization language, the field ip represents an operator 

^(x,t)-^a„^„(x)e-^" t , (5) 

n 

where the ip n (x) are eigenstates of the Schrodinger equation, the u) n are the 
corresponding eigen-energies, and the a n are the corresponding annihilation 
operators for particles in that mode. If we specialize to the case where there 
is no external potential [V(x) = 0], then this becomes 

$(x,t) -» J a k e lk ' x -^ fc *, (6) 

which looks just like the quantization of field in terms of plane waves that 
you're used to from relativistic quantum field theory. 

Recall that in single-particle QM, ip*ip gives you the probability density. 
In second-quantized QM, the analogous statement is that gives you the 
number density, so 

N= f ^ = ^>ta n . (7) 

To describe a system of particles with a given number density n, it's convenient 
to use the grand-canonical ensemble and introduce a chemical potential term 
fiN in the Hamiltonian or Lagrangian. So, our final Lagrangian for a free 
non-rclativistic Bose gas is 

L = J d 3 x V>* (idt + V 2 + ,u — V(x)^j i>. (8) 
4 Interactions 

Now let's include a two-body potential U(x — y) between the atoms (which 
is adequate in the dilute limit). Remembering that tp*ip(x) is the number 
density at x, we need 

imt = -\\ V*V>(x)t/(x-y)V*V(y). (9) 



If U is localized, then, in the low-energy limit (wavelength 3> a), we can 
replace it by an effective ^-function: 



tf(x-y)-— <J(x-y), 
m 



(10) 



where a is the scattering length and can be measured for a given type of 
atom. For this substitution to be valid, the typical inter-particle spacing n -1 / 3 
(where n is the number density) should be large compared to the scattering 
length. In leading non-trivial order in this "diluteness" expansion a <C n" 1 / 3 , 
the Lagrangian for the interacting Bose gas is then simply 

V>* (idt + 7^7 V 2 + n - V(x)) 4, - —(n>) 2 
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(11) 



I will henceforth focus on the case of zero external potential V. In that 
case, this looks just like relativistic 4> A theory, but with c/)*(~df)(f> replaced 
by ip*idtip. In particular, — fj, now plays the role that m 2 ff (T) did in the 
introduction, and one can see that there will be a phase transition depending 
on whether — /i(T) > (symmetry restored phase) or — /i(T) < (symmetry 
broken phase). In the latter case, the expectation (ip) ^ is what's known as 
the Bose-Einstein condensate. Q 

5 Now Do the Usual 



To study the phase transition, go to Euclidean-time formalism. As in the in- 
troduction, note that all but the zero-Euclidcan-frequency modes ipo decouple 
at large distance, leaving 



dr L\ 







d 3 a 



V-'o 



~ V 2 - Meff + V{x) ] ^0 + — (^o) 2 
Am / m 



(12) 

[Corrections to this effective theory turn out to be higher order in the dilute- 
ness expansion.] This Lagrangian can be made to have standard field-theory 
normalization for the kinetic term by rescaling ipo — §\j2mT: 



dr L 



E 



(i 



(3 / d 3 x [|V0| 2 +7#| 2 +^| 4 ] , 



(13) 



where r = —2mfj, e g and u = AnamT. Note that the diluteness condition 
a <C n 1 / 3 implies that the coupling u should be considered as "small," since 



a In graduate statistical mechanics, you learned for free particles that the condensate is 
the particles in the k = mode. The connection is that the contribution of k = to the 
number density n = (tp*ip) is no = (tp)*(ip). 



it is proportional to a. So, we have here an 0(2) field theory in 3-dimensions 
(since (j> is complex) which is weakly-coupled at short distances]^ 

6 A Goal: Calculate T c as a Function of n 

An interesting thing to try with this effective theory is to calculate the correc- 
tion, due to interactions, to the ideal gas result for the Bose-Einstein conden- 
sation temperature T c . Actually, it turns out to be technially slightly easier 
to calculate the shift An(T) in the critical density due to interactions, at 
fixed temperature, rather than the shift AT(n) in the critical temperature, 
at fixed density. The two are easily related. Then recall that the density 
n = ~ (0*0)- If one were to do perturbation theory, the sort of dia- 

grams one would calculate for An would then be of the form 




where the cross represent the operator 4>*4>. But perturbation theory breaks 
down at the transition; so what to do? n 
One possible technique, implemented by Baym, Blaizot, and Zinn- Justing 
is to try the large N approximation for solving the 3-dimensional theory, 
setting N—2 at the end. At leading order in large N, the graphs which 
contribute to (0*0) — ► ((ft ■ 4>) are 




Baym et al. find T c = Tq(L + 2.33 an 1 / 3 ) plus higher orders in 1/N and in 
an 1 / 3 . Boris Tomasik and 10 have analyzed the next-order corrections in 1/N 
and find that they change the coefficient 2.33 by only -26% for N=2. This 
correction is surprisingly small and suggests that large N might not be too 
bad for T c \ 
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b More technically, this means weakly coupled at the short-distance scale k ~ \J1raT where 
it breaks down, where one wants to match it to the original theory. 
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